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Abstract 

The motivation of this paper is to introduce a novel framework for the restoration of images corrupted 
by mixed Gaussian-impulse noise. To this aim, first, an adaptive curvelet thresholding criterion is proposed 
which tries to adaptively remove the perturbations appeared during denoising process. Then, a new statistical 
regularization term, called joint adaptive statistical prior (JASP), is established which enforces both the local 
and nonlocal statistical consistencies, simultaneously, in a unified manner. Furthermore, a novel technique 
for mixed Gaussian plus impulse noise removal using JASP in a variational scheme is developed—we refer 
to it as De-JASP. To efficiently solve the above variational scheme, an efficient alternating minimization 
algorithm is developed based on split Bregman iterative framework. Extensive experimental results manifest 
the effectiveness of the proposed method comparing with the current state-of-the-art methods in mixed 
Gaussian-impulse noise removal. 


Index Terms 

Mixed Gaussian-impulse noise removal. Adaptive curvelet thresholding, nonlocal self-similarity. 

I. Introduction 

Image denoising as a fundamental problem in image restoration has been extensively studied in image 
processing and computer vision (see [1]-[12] for few various representative works). In this paper, we 
address the problem of image restoration from mixed Gaussian-impulse noise data. Suppose the clean 
image u e is corrupted by both additive white Gaussian noise (AWGN) of p e (with variance 
cr^) and an impulse noise. Also, let Nimp represent the process of image degradation with impulse noise. 
Accordingly, the overall degradation by mixed Gaussian-impulse noise can be modeled as the following: 

/ ~ V)"! (D 

where / G is the observed image (noisy measurements). In the literature, there are two common types 
of impulse noise: salt-and-pepper (SP) noise and random-valued (RV) noise. Let denote Ua^b as the pixel 
value of image u at location (a, b) G X, and [timin, t/max] be the dynamic range of u. The operator Nimp{') 
is defined as follows: 

{ da^bi with probability r/2, 

with probability r/2, (2) 

with probability 1 — r, 
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where da^b and ^ are in the range of [dmin, <^max] and independent from Ua^b- When da^b = dmin and 
d'ab — <^max, the model is SP noise; however, if da,6 = ^ is identically and uniformly distributed 

random numbers in [dmin, c^max] with probability r, the model is called RV noise. Median-type filters, 
e.g. adaptive median filter (AMF) [13] and adaptive center-weighted median filter (ACWMF) [14], are 
widely used in the literature to remove impulse noise. 

Let AT denote the set of locations of noisy pixel candidates corrupted by impulse noise, and accordingly, 
the set of pixels that are likely to be uncorrupted by impulse noise is defined SiS A = X/Af. Therefore, 
the degradation matrix $ can be defined as ^a,6 := {1 : (a, 6) G Al or 0 : (a, 6) G JV} Q More precisely, 
the pixels probably being corrupted by impulse noise are assigned with 0 while the rest of pixel values 
are more likely to be corrupted by AWGN. Delving further into this issue reveals this fact that, ^ o / 
follows an approximate Gaussian distribution, where o denotes the element-wise multiplication between 
two matrices. Also, due to the ill-posed nature of (1), regularization-based methods are often used to 
convert (1) into a well-posed problem by imposing a priori information of the underlying signal [7], 
[9]-[12]. Hence, from the Bayesian inference, the following regularization-based framework for mixed 
Gaussian-impulse noise removal is written as: 

T^A\\^o{f-u)\\%+XJ{u), (3) 

where the first term is a penalty that represents the closeness of the solution to the observed scene; the 
second term is a regularization term which represents a priori information of the original scene; and A is 
a regularization parameter that balances the contribution of both terms. The term J{u) could be various 
choices, e.g., Tikhonov regularization [15], Geman and Reynolds’ half quadratic variational model [16], 
Rudin, Osher and Fatemi’s total variation model [17], Mumford-Shah model [18], and framelet based 
model [7]. These approaches do not need to find the damaged pixels and perform well in impulse noise 
removal. However, for the case of images corrupted by mixed Gaussian impulse noise, the Gaussian 
noise is not treated properly [12]. In recent works, the local smoothness and the nonlocal self-similarity 
properties (or image sparsity prior) exhibited in natural images are combined into the final cost functional 
of image restoration solution to achieve better performance, e.g., see [10]-[12]. 

The main contributions of this paper are listed as follows. First, we propose an adaptive curvelet 
thresholding (ACT) criterion, in a statistical manner, for adaptively characterizing the perturbations ap¬ 
peared during noise removal process—exploiting local consistency. Second, we establish a new statistical 
regularization term, called joint adaptive statistical prior (JASP) which enforces both the local and nonlocal 
statistical consistencies, simultaneously, in a unified manner. Third, we propose a novel technique for mixed 
Gaussian plus impulse noise removal using JASP in a variational scheme, we refer to it as De-JASP. To 
solve the above variational scheme, an efficient alternating minimization algorithm is developed based on 
split Bregman iterative framework [19]. 

The rest of this paper is organized as follows. Section II provides a brief background on a recently 
proposed modeling for nonlocal self-similarity. The proposed ACT is introduced in Section III. The 
proposed regularization term, along with its incorporation into the variational framework of mixed noise 
removal, and implementation details are introduced in Section IV. Numerical results are given in Section 
V and finally. Section VI concludes the paper. 

H. Background on Nonlocal Statistical Modeling (NLSM) 

Motivated by the success of BM3D [4] in image denoising, Zhang et al. [11] proposed a model for 
nonlocal self-similarity prior information, which is employed efficiently in image restoration. The NLSM 

Tn this paper, similar to [5], [7]-[12] and for the sake of a fair comparison, $ is determined by AMF and ACWMF for SP and 
RV noise detection, respectively. 
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Fig. 1: (a)-(b) The histograms of two successive orientations at one scale of the curvelet coefficients of the clean image “Fence” 
(which is shown in Fig. 2(a)). (c)-(d) The histograms of the same orientations at the same scale of the curvelet coefficients of the 
initiated image uinit (which is shown in Fig. 2(c)), where, it is reconstructed by employing AMF on noisy image of Fig. 2(h). 
(e)-(f) PDFs of noise distribution (empirical perturbations) in natural log domain at the same orientations and scale (corresponding 
to the histograms of Figs. l(a)-l(d)) along with their best fitted PDFs: the noise distribution (blue), the Gaussian fit (red), and the 
Laplacian fit (black). 


explores the nonlocal self-similarity by means of the distribution of the transform coefficients, which are 
obtained by transforming the 3D array generated by stacking similar image patches. To elucidate on, First, 
the image u (of size s/n x ^/n) is divided into P overlapped patches of equal sizes, i.e., Up^, (of size 
^/B^ X ^/Bp) at location where /' = 1, 2, • • • , P. Then, for each patch the best c matched patches are 
found within a searching window (of size cc x cc). Next, these similar patches are stacked into a 3D array, 
which is called a group. Now, a 3D transform, T^^('), is conducted on the group to obtain its 
transform coefficients, followed by arranging them in a lexicographic order. The mathematical formulation 
of the NLSM for self-similarity in transform domain can be written as: 

p 

i'=i 

where ^nlsm(') corresponds to the NLSM operator. To put it simply and briefly, after obtaining 0^, the 
new estimate of u is achieved by u = ^nlsm(®u )5 where f^NLSM(') corresponds the inverse operator of 

^ NLSM • 


III. The Proposed Adaptive Curvelet Thresholding (ACT) 

The curvelet transform [20] is a multi-scale directional transform which provides a new multi-resolution 
representation with several features that are superior to existing representation, i.e., wavelets and steerable 
pyramids, thereby gaining good performance in image denoising [1], [21]. For these reasons, here, an 
adaptive denoising/shrinkage procedure in curvelet domain via estimation of clean curvelet coefficients 
from available rough estimated data is proposed 

We assume that the initially estimated image Uinu achieved by the median filter (AMF or AC WMF) 
is perturbed and model that perturbation by additive noise, with no prejudging on being AWGN or any 
specific species of noise. To elucidate on, this modeling of perturbations by additive noise helps us to model 
the noise probability density function (PDF) effectively and being able to solve the ensuing optimization 
problem efficiently. Suppose, the initiated image coefficients in curvelet domain ('d = can be 

^In matrix notation, we can write d = and u = where d is the curvelet coefficient of function n, is a forward 

curvelet transform and is its Moore-Penrose pseudo-inverse transform (corresponding to the minimal dual 

synthesis frame). Due to the Parseval tight frame property of curvelet transform (i.e., = I), we have 
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written as: 

+ ^i,z,k, (5) 

where ^ and ^ are the initiated image curvelet coefficient (available rough coefficient which is assumed to 
be noisy) and the noise-free (or clean) curvelet coefficients, respectively, and n is unknown additive noise 
(or perturbation) vector. Also, the three indices of j,/,k are curvelet triple index of scale, rotation and 
position, respectively. The problem of image denoising in curvelet domain can be expressed as estimation 
of clean coefficients from noisy data with Bayesian estimation techniques. Either minimum mean-squared 
error (MMSE) or maximum a posteriori (MAP) estimator is used for solving this problem, the solution 
requires a prior knowledge about the distribution of curvelet coefficients. More specifically, these methods 
are optimized w.r.t. the marginal statistics of the coefficients within each sub-band, by imposing a prior 
distribution on the clean transform coefficients [6]. 

Eigs. 1(a)-1(b) show the histograms of two successive orientations at one scale of the curvelet coefficients 
of the clean image “Eence” (256 x 256)—see Eig. 2(a). Obviously, it can be observed that the distributions 
are characterized by a very sharp peak at zero amplitude and extended tails to both sides of the peak 
(leptokurtic distribution). This leptokurtic property is observed in all histograms of each scale and orien¬ 
tation of all test images. In this paper, Laplacian distribution is chosen to model the marginal distributions 
of curvelet coefficients of the clean image, i.e., Pr^(^j^^^k) = (l/>/2crj,^,k) • exp ( |'^j,Z,k|/c^j,Z,k)5 

where / k is the standard deviation of clean coefficients. This choice is despite the fact that marginal 
distributions of curvelet coefficients have significantly heavier tails than a Laplacian distribution; tend to 
go zero slower than a Gaussian distribution; and can be well modelled by a hyper-Laplacian distribution. 
But, actually, this choice makes a trade-off between modeling the coefficients statistics accurately and 
being able to solve the ensuing optimization problem efficiently. 

The histograms of the same two successive orientations at the same scale of the coefficients of the 
initially estimated image Uinu —see Eig. 2 (c)—are shown in Eigs. 1(c)-1(d). More specifically, these 
histograms are obtained using the curvelet coefficients of the initiated image Uinit, which is achieved by 
the initial median filtering step. As can be seen, here, the leptokurtic behavior and heavy tails do not exist 
in the shown distributions. More accurately, these two properties do not exist in all distributions of each 
scale and orientation of all initiated images achieved by this step. 

Considering (5), the noise distribution can be obtained using the curvelet coefficients of the clean image 
and the initiated image. Eigs. 1(e)-1(f) show the PDEs of noise distribution at the same two successive 
orientations of the same scale corresponding to those of Eigs. 1(a)-1(d). It can be observed that the 
empirical distributions of noise coefficients can be well characterized by Gaussian distributions, 

while the Laplacian distributions have much larger fitting errors. This observation motivates us to model 
the PDE of nj^/^k by a standard Gaussian distribution: Pr^(nj^/^k) = {^/V^cFuj i) - exp k/2<^nj J’ 

where ^ is the noise variance in curvelet domain, which can be estimated using robust median estimator 
[22] and Monte-Carlo simulations [20]. 

Consider the marginal distribution of the clean curvelet coefficients. Since in image denoising problems, 
we do not have access to the clean image, the variance of the clean image curvelet coefficients ^ k can 
be estimated by maximum likelihood (ML) estimator, through averaging over the neighboring initiated 
image coefficients (available ^j,^,k) within a squared window [23]. 

In order to estimate ^j,z,k from the available observation the MAP estimator is employed: 

= argmax \ Pr 



( 6 ) 
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(a) (b) (c) (d) 


Fig. 2: (a) Original Fence image (of szie 256 x 256). (b) The noisy image which is corrupted by Gaussian plus salt-and-pepper 
impulse noise with a = 30 and Csp = 30% (PSNR=9.92 dB, SSIM=0.143). (c) The initiated image using AMF (PSNR=18.81 dB, 
SSIM=0.432). (d) The recovered image via ACT by 11 iterations (PSNR=22.11 dB, SSIM=0.589). 


By considering (5) and using Bayes’ rule, we get: 


= arg max 

'’^3,IF 


= arg max 

^3,IF 


Pr„ (i?i,i,k - i?i,i.k) • Pr^ (^?i,i.k) j- 


(7) 


Therefore, (7) allows us to write this estimation in terms of the PDF of the noise coefficient (Pr^), and 
the PDF of the clean image coefficient (Pr^). By substituting Pr^ and Pr^ into (7) and doing some 
manipulations, the estimation of is achieved as: 


= sign(i?j,;,k) max (O, |^j,i,k| 




( 8 ) 


where, in fact, (8) is the classical soft shrinkage function [24] defined as Sp{y) := sgn(^) •max(0, \y \ — p). 
Not only can the proposed ACT applied for the initiated image achieved by the initial median filtering 
step, but also it can be employed in further iterations—by using the last approximated image achieved 
in previous iteration as the noisy image for the current iteration—which leads to a simple and effective 
mixed Gaussian-impulse noise removal scheme, e.g., see Fig. 2(d). 


IV. The Proposed Image Denoising via JASP (De-JASP) 

The proposed joint adaptive statistical prior (JASP) is defined by integrating both the local information 
prior (which depicts the local smoothness and geometric regularity of image structures achieved by discrete 
curvelet transform (DCuT) [20]) and the nonlocal self-similarity prior (corresponding to the NLSM in 3D 
transform domain achieved by the method introduced in Section II). In a mathematical expression, the 
proposed JASP is expressed as: 

'^JASp(^) ^DCuT (^) "^^NLSM (^) 

= ||^«|k+r||e„|k, (9) 

where the first and the second terms represent the image local smoothness prior and nonlocal self¬ 
similarity prior, respectively. Also, r is a regularization parameter which controls the trade-off between 
two competing (statistical) terms. 
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By substituting the proposed JASP (9) for the regularization term J{u) in the regularization-based 
framework of (3), and by introducing two auxiliary variables ^ and w, the proposed optimization problem 
for image recovery is expressed as follow: 

mm{h^oif-u)\\l+X\\ne,+r'\\eM 

u,ij,w Z 

S.t. = U = W, W = (10) 


In order to solve the above minimization problem, an alternating split Bregman iterative algorithm [19] 
is invoked. We finally obtain the following schemes: 

(y/=+i^^fe+i^y,/=+i) min |^||^ o {/- u)||2 + A||t?||<!j 

u,'U,w k Z 

+r'||e»lk. +^\\^-'iu- b’^Wl + l^\\u-w- 0 %^} (11a) 

j^k+l _ ^k+1 ^^fe+1 (1 lb) 

cfc+l ^^/=_y/=+l^y,/=+l_ (ll^.) 


Here, /ii and 112 are fixed value parameters for improving the numerical stability of the algorithm. 

Given and the u sub-problem of (11a) consists of minimizing a strictly convex quadratic 
function that can be solved easily, i.e., ^ where p = ^ p 2 and V = 

— 6^) + /i2(^^ + c^) + o /). In order to avoid computing the matrix inversion in the first 
term, the matrix inversion lemma of the Sherman-Morrison formula is applied to that term. Thus, the u 
sub-problem can be obtained as: 


yfe+i ^ 


+ ° ■ 


( 12 ) 


By given u in hand (according to (12)) and considering ^ + b (for simplicity, the superscript k 

is dropped without confusion), the ^ sub-problem of (11a) becomes: 

(13) 

V Z PjI 

By these transformations, we regard as some type of the noisy unknown coefficient More precisely, 
the sub-problem of (13) can be interpreted as the denoising in curvelet coefficient domain. Considering the 
fact that the unknown variable ^ is component-wise separable in curvelet domain, each of its component 
can be obtained by a component-wise shrinkage procedure: 


= arg min | 1, 

( z p^l ) 


(14) 


where its solution is the same as the one introduced in (8) with regarding To elucidate 

on, the ^ sub-problem of (13) can be interpreted as the proposed ACT introduced in Section III. Thus, 
the parameter A is self-adaptive. 

Given u, analogously, the w sub-problem of (11a) can be written as: 

min - 2;||| +—||0„ 11^1, (15) 

w Z P 2 

where z = u — c (for simplicity the superscript k is omitted without confusion). Owing to the complicated 
definition of it seems difficult to solve (15) directly. In order to solve (15) amenable, according to 
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TABLE I: The Proposed De-JASP Framework 


Input: /, Bp.uj.r', 111, ii2,Tol, maximum iteration number 

Output: n*: Recovered image; 

Initialization: Set A: = 0; (6°, c°) = (0, 0) 

Uinit = median-filtering(/); compute 

While a stop criterion is not satisfied do 

1 . = u = u^\ + b^\ 

2. solve -i? sub-problem using (14) and concatenating all coefficients 

3. solve w sub-problem using (19) 

4. solve u sub-problem using (12) 

5. compute SSIM(n^+^, u) 

6. compute dijf = abs(s^+^ — s^) 

7. update using (11b) and (11c), respectively 

8. k i — A: + 1 
end While 

stopping criterion: k = k'^^^ or diff < Tol. 


[11], a reasonable assumption is used which leads to obtain a closed-form solution of (15). By assuming 
that all elements of u; — 2 : (such that w^z e M^) are independent and identically distributed (i.i.d.) with 
zero-mean and variance and by invoking the law of large numbers in probability theory, for any 

e' > 0, it leads to: 

lim Pr{|-||«; - z\\j - (t'^| < ^} = 1. (16) 

n^oo ' n Z 

Due to the orthogonal property of 3D transform and denoting the error vector in 3D transform 

domain by ©e = ©^ — Qz (such that Oz G where Kq = Bp * c * P), it is concluded that all 

elements of ©e are i.i.d. with zero-mean and variance . Thanks to the law of large numbers, by doing 
the same manipulations with (16) to ©g, for any e' > 0, it yields: 

Therefore, according to (16) and (17) the following relationship is described (almost surely): ^11^7; — 41 = 
■^W^w — ©ZIII 2 ’ where incorporating it into (15) leads to: 

argminl||0^ - ||0„,. (18) 

Z /7'2 77/ 

Since the unknown variable ©^^ is component-wise separable in (18), each of its component Qw{j')^j' = 
1, • • • ,P© can be independently obtained by a component-wise (soft) shrinkage procedure. Thus, the 
closed-form solution of (15) is written as: 

^ ~ ^nlsm(Qw) ~ (Q;g))- (1^) 

The proposed De-JASP algorithm is delineated in Table I. 

^It is worth emphasizing that the above assumption does not need to be Gaussian, Laplacian or generalized Gaussian process. 








TABLE II: PSNR/SSIM COMPARISONS FOR Gaussian plus Salt-and-Pepper Noise Removal. 



Algorithm 

Test Image 

r* sp 

cr 

Noisy Image 

Init. Image 

IFASDA [7] 

JSM [11] 

De-JASP 


0.2 

20 

12.10/0.147 

22.48/0.529 

24.92/0.756 

27.73/0.795 

28 . 06 / 0.835 

Barbara 

0.2 

30 

11.69/0.133 

19.57/0.385 

20.46/0.514 

22.76/0.537 

24 . 83 / 0.676 


0.3 

30 

10.21/0.091 

19.10/0.364 

20.34/0.502 

22.63/0.533 

24 . 43 / 0.655 


0.2 

20 

12.14/0.134 

23.08/0.510 

27.81/0.771 

28.40/0.784 

29 . 22 / 0.833 

Boat 

0.2 

30 

11.72/0.122 

19.93/0.373 

22.74/0.506 

23.22/0.520 

25 . 88 / 0.679 


0.3 

30 

10.24/0.085 

19.47/0.355 

22.67/0.499 

23.15/0.518 

25 . 47 / 0.659 


0.2 

20 

11.77/0.212 

21.58/0.560 

24.99/0.782 

26.60/0.791 

26 . 79 / 0.812 

Fence 

0.2 

30 

11.39/0.197 

19.13/0.443 

21.03/0.577 

22.35/0.583 

24 . 15 / 0.693 


0.3 

30 

9.92/0.143 

18.81/0.432 

20.82/0.569 

22.16/0.577 

23 . 80 / 0.679 


0.2 

20 

11.89/0.095 

23.15/0.407 

27.78/0.759 

28.83/0.765 

29 . 65 / 0.837 

Parrot 

0.2 

30 

11.50/0.087 

19.97/0.277 

22.65/0.424 

23.51/0.438 

26 . 21 / 0.640 


0.3 

30 

10.02/0.063 

19.51/0.264 

22.46/0.411 

23.31/0.435 

25 . 61 / 0.616 


V. Experimental Results 

In this section, we evaluate the performance of the proposed method—De-JASR The performances 
of our experiments are evaluated on 4 gray-scale test images of Barbara, Boat, Fence, and Parrot 
(of size 256 x 256). To evaluate our simulation results, we use two applicable quality assessors, the 
pick signal-to-noise ratio (PSNR) in dB and the structural similarity (SSIM) [25]. All the experimental 
results reported in this paper are averaged over 10 independent trials. Also, the second generation of 
DCuT via wrapping [20] is employed in our implementations. To demonstrate the effectiveness of the 
proposed method in image denoising, we have compared it with two competitive mixed noise removal 
techniques: (i) iterative framelet-based approximation/sparsity deblurring algorithm IFASDA [7], and (ii) 
image denoising via joint statistical modeling (JSM) [11]. It is worth emphasizing that these two recently 
proposed denoising methods are among the state-of-the-art techniques in mixed-noise removal which 
achieve the best performance so far. Note that, to make a fair comparison among the competing methods, 
we have carefully tuned their parameters to achieve the best performance. Also, for the sake of fair 
comparison, the same test conditions are used in all experiments; i.e., the same noisy measurements are 
applied for each method. Due to the space limitations, only parts of the experimental results are shown 
in this paper. Interested readers may contact the corresponding author for all the images and the source 
code. 

In the following, the effectiveness of De-JASP is investigated. In our implementation, in process of 
finding NLSM for self-similarity in 3D transform domain, the size of each patch is set to 8 x 8. The 
size of training window for searching matched blocks is set to 20 x 20, and the number of best matched 
blocks is set to 10. The other main parameters of the proposed algorithm are empirically set: — ^0, 

Tol = 10“^, til = 0.16//, //2 = 0.84/i and t^ = 2.7 x 10“^. Also, the orthogonal 3D transform denoted 
by is composed of 2D DCT and ID Haar transform. 

Table II lists the PSNR and SSIM values of the recovered noisy test images obtained by the proposed 
De-JASP compared to those obtained by JSM and IFASDA in different noise level hypotheses (best results 
are emphasized in bold face). From Table II, it can be inferred that the performance of the proposed De- 
JASP is superior to those of compared algorithms in terms of both employed objective quality assessors. 
The average PSNR (SSIM) gain of De-JASP over JSM and IFASDA, in mixed AWGN-\-RV noise removal, 
can be up to 2.31 dB (0.154) and 3.55 dB (0.167), respectively. For visual comparison, some recovered 
images by the competing methods, in different noise levels, are shown in Figs. 4-7. It can be observed 
that De-JASP outperforms the other competing methods and reproducing clearer images. Better visual 
comparison can be made by zooming the images on the screen. 






Fig. 4: Visual quality comparison of mixed Gaussian-impulse noise removal on image “Barbara”, (a) noisy image corrupted by 
AWGN+SP noise with a = 20 and rsp = 20%. Denoised results by (b) IFASDA [7]; (c) JSM [11]; and (d) the proposed De-JASP. 
For numerical comparison, see Table II. 


Fig. 5: Visual quality comparison of mixed Gaussian-impulse noise removal on image “Boat”, (a) noisy image corrupted by 
AWGN+SP noise with a = 10 and rsp = 60%. Denoised results by (b) IFASDA [7] (PSNR=27.58 dB, SSIM=0.769); (c) 
JSM [II] (PSNR=28.06 dB, 0.793); and (d) the proposed De-JASP (PSNR=28.37 dB, SSIM=0.8I4). 


(m) (n) (o) (p) 

Fig. 6: Visual quality comparison of mixed Gaussian-impulse noise removal on image “Parrot”, (a) noisy image corrupted by 
AWGN+RV noise with cr = 20 and Vrv = 20%. Denoised results by (b) IFASDA [7] (PSNR=26.34 dB, SSIM=0.726); (c) JSM 
[II] (PSNR=26.96 dB, 0.789); and (d) the proposed De-JASP (PSNR=27.42 dB, SSIM=0.8I6) 
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(q) (r) (s) (t) 

Fig. 7: Visual quality comparison of mixed Gaussian-impulse noise removal on image “Fence”, (a) noisy image corrupted by 
AWGN+SP+RV noise with a = 10, Vsp = 10% and Vrv = 10%. Denoised results by (b) IFASDA [7] (PSNR=20.11 dB, 
SSIM=0.516); (c) JSM [11] (PSNR=21.79 dB, 0.622); and (d) the proposed De-JASP (PSNR=23.19 dB, SSIM=0.698) 


VI. Conclusion 

In this paper, first, an adaptive curvelet thresholding (ACT) criterion is introduced to adaptively char¬ 
acterize and abate the perturbations which are appeared during denoising process. Then, a novel joint 
adaptive statistical prior (JASP) and a new strategy for mixed Gaussian-impulse noise removal via JASP, 
called De-JASP, is proposed. Extensive experimental results clearly confirm that De-JASP significantly 
outperforms the current state-of-the-art mixed Gaussian-impulse noise removal techniques. 

For future work, we would like to propose an effective nonlocal self-similarity modeling, and pursue 
this direction to extend our proposed statistical prior to other image restoration applications such as image 
deblurring and inpainting. 
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